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ABSTRACT 


Adhesive bonding has been used to join or repair metallic and composite structural 
components to achieve or restore their designated structural stiffness and strengths. However, 
current analysis methods and empirical databases for composite bonded patch repair or joints 
are limited to flat structures, and there exists a very limited knowledge on the effect of 
curvature on the performance and durability of composite bonded joints and repairs. Recently, 
a novel finite element formulation was presented for developing adhesive elements for 
conducting 2.5-D simplified stress analysis of bonded repairs to curved structures. This report 
presents the work on optimal design of bonded composite patch using the newly developed 
adhesive element. Three types of design variables considered separately are (a) geometrical 
boundary shape of bonded composite patch, (b) ply drop-offs contour of bonded composite 
patch, and (c) fiber orientation of each ply in bonded composite patch. Sequential Linear 
Programming (SLP) method and one-dimensional direct search method are employed as the 
optimization algorithm in conjunction with a fully implemented mesh generation algorithm 
into which new features were incorporated. Several different objective functions were 
proposed to optimize the bonded patch. 
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1. INTRODUCTION 


Adhesive-bonded patching is one of the widely used repair techniques to cracked or damaged 
metallic and composite structures [1,2]. In this technique, a composite patch is bonded to the 
parent structure to reinforce the cracked zone [3] with the intent of restoring the structure to 
its original design specifications. This technique has successfully addressed some of the 
aging aircraft problems [1], and there exists a large number of research publications available 
on analysis, testing and design of adhesive-bonded patching repairs. A bonded patch is 
similar to an adhesive-bonded joint, thus the concepts and methodologies developed for 
analyzing stresses in bonded joints can be readily applied to conducting stress calculation for 
a bonded patch. There is a large amount of information on stress analysis of the bonded joints 
and repairs in the literature. Most of the currently available analysis methods and empirical 
databases for composite bonded joint design and for composite patch repair designs are 
limited to flat plate and/or flat laminate geometries [1, 4-6]. Compared to the knowledge of 
bonded joints/repairs to flat surfaces, current knowledge on the effect of curvature(s) on the 
performance and durability of composite bonded joints and bonded composite repairs is 
limited. In the first phase of this project, the effect of curvature on both peel and shear 
stresses in adhesive layer in adhesive-bonded curved patches was investigated by using three 
newly-developed adhesive elements [7,8]. In the finite element formulations, it was assumed 
that the peel and shear strains be constant across adhesive thickness. Some illustrative 
examples were also presented to demonstrate the effect of curvature on peak peel and shear 
stresses in the adhesive layer under selected loading and boundary conditions. 

The objective of this phase of the project is to develop optimization scheme for optimum 
design of bonded patches. Considered are three separate sets of design variables, namely, the 
geometrical shape of a patch, ply drop-off (tapering) scheme and fiber orientations for each 
ply. Sequential Linear Programming method and one dimensional direct search method are 
used to optimize independently the three sets of design variables. The objective function is 
selected as the maximum peel stress, maximum shear stress, and maximum Von Mises stress, 
respectively. To facilitate mesh changes in the optimization process, an automatic mesh 
generator is developed to achieve a minimum the influence of shape changes on the peak 
objective stress. To take into account the ply drop offs, a new adhesive element scheme is 
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introduced to allow variation of neutral plane from element to element, which was validated 
by comparing with a full 3-dimensional finite element analysis using a commercial code. 

2. STRESS ANALYSIS USING ADHESIVE ELEMENTS 

Figure 2.1 depicts a typical curved surface with a bonded patch from one side. To optimize 
the bonded patch, it is essential to employ a simple and effective finite element analysis 
procedure. In this report, the 2.5D adhesive elements developed by the authors [7,8] are 
employed. In this section, a brief review the employed elements is given and a discussion on 
formulations for calculating nodal stresses in the adhesive layer is also presented. 



Figure 2-1 A curved surface with a bonded patch 

2.1. Review of degenerated shell elements 

There are two kinds of shell element for structural analysis, namely, plate-shell element and 
degenerated shell element. For the former one, i is based on plate element combining a 
drilling degree to simulate free rotation in 3-D space. However, plate-shell element is limited 
to cases in which all the nodes in an element must be located in one single flat plane. This 
introduces extra extreme difficulty to automatic mesh generator for those elements with more 
than 4 nodes. Although triangle plate-shell element is more attractive for this circumstance, it 
needs much more fine mesh to gain a satisfactory precision of stress analysis. For degenerated 
shell element, every node has its own nodal coordinate system and it just needs to put all the 
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nodes of an element in one single curved plane, which does not add any limitation to auto¬ 
meshing. Therefore, 8-node degenerated shell element with the first order shear deformation 
is selected to analyze both curved surface and bonded patch. 

According to the first order shear deformation theory, Ihe nodal displacements of a shell 
element can be defined in term of the displacements in its mid-plane in the local coordinate 
systems 

U(x, y,z) =u°(x,y)+ zQ y (x,y) 

V(x, y, z) = v° (x, y) - z.Q x (x, y) (2-1) 

W(x, y,z) = w(x,y ) 

where w°, v°, w are the translational displacements in the mid-plane of the shell, and Q x , 0 are 
the rotations of the directional normal about x and y coordinate respectively. In the shell 
element, each node has 5 degrees of freedom, i.e., 3 translational degrees of freedom u°, v°, w 
and 2 rotational degrees of freedom 0 t , 0 v . 

2.2. Review of adhesive element 

Figure 2-2 depicts a 16-node adhesive element, which consists of one pseudo-brick element 
modeling the adhesive layer and connecting the two shell elements. For the pseudo-brick 
element modeling the adhesive layer between the two shell elements, the following three 
basic assumptions were employed according to the characteristic of the adhesive layer: 

• The patch is perfectly bonded to the curved surface with uniform thickness and no 
debonding occurs in the whole adhesive layer. 

• The thickness of adhesive layer is very thin comparing to that of patch and curved 
surface. Then the strains in the adhesive layer are assumed to be uniform across the 
thickness of the adhesive layer. 

• Only three out-of-plane stresses a z ,J iz , J xz are considered and other three in-plane 
stresses CT a , ct y , y are assumed to be equal to zero. 
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Figure 2-2 A 16-node adhesive element 

Based on the basic assumptions above, the peel and two shear strains in the adhesive between 
pseudo-node i' and j' (see Figure 2-3) are constant and can be expressed as follows: 


< =-(w, ~ w ,) 


1 


h 


y; = - (Y- - v r )=- (v, u - v; + -^e,+ -^e xj ) 

■ u° — . — — 0 .) 

1 2 ' 2 » 


Y.=7(^-«/) = 7K° 


( 2 - 2 ) 


where the superscript “°”denotes the mid-plane, subscripts i and j denote the upper and lower 
shell elements, respectively. h u and hi represent the thickness of upper and lower shell 
elements, t is the thickness of the adhesive layer. 
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2.3. Stress interpolation 


For an iso-parametric displacement element, the stresses at the Gaussian points have better 
precision, while those on the node don’t posses good precision because displacements and 
stresses in an element have different order of continuity. Then the nodal stresses can’t be 
obtained directly with a good accuracy. For bonded patch/joints, there always is a stress 
concentration near the boundary of the adhesive layer. This indicates that the stresses along 
the boundary are the key parameters to determine bonding strength of the adhesive layer. To 
overcome this drawback, a stress interpolation scheme is employed to gain the nodal stresses 
based on the stresses at the Gaussian points [9]. 


a 

b 

c 

b 

■s/3 


b c b 
a b c 
b a b 
c b a 


'in 


1 , V , 1 s/3 

1H-, b = —, c = 1- 

2 2 


(2-3) 


where a,, a 2 , a 3 , a 4 are the interpolated nodal stresses, ct ( , a n , <j m , a IV are the stresses on 
the Gaussian point. The locations of node and Gaussian point are schematically shown in 
Figure 2-4. 



Figure 2-4 Locations of node and Gaussian point in an 8-node element 
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3. SHAPE OPTIMIZATION OF A BONDED PATCH 


3.1. Shape optimization method 

Sequential Linear Programming (SLP) is an effective method for shape optimization problem 
[10]. It concerns on linearising the objective function and constraints by sensitivity analysis 
and then the nonlinear optimization problem being solved by the linear programming method. 
In this section, we give a brief description of SLP first and then present the detailed 
formulations for shape optimization of the bonded patch. 

3.1.1. Sequential linear programming 

A structural shape optimization problem can be described as the minimization of an objective 
function subject to a number of prescribed constraints. It may be stated mathematically as 


Minimize f(X) 
Subject to 


| h i (X) = e i (i = l,...,p) 
\gj(X)<b j (j =!,...,/«) 


(3-1) 


Where f(X) is the objective function, h t {X) and g j(X ) represent the equality and inequality 

constraint functions, X is the design variable vector which contains n design variables, and p 
and m are the total number of equality and inequality constraints, respectively. 

The optimization formulation as above is nonlinear in general. The linearization of the 
formulation is a must if the linear programming is going to be employed to solve the problem. 
The linearized formulation becomes 


a a/ 

Minimize ) -—Ar 


k =i 


k 
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Subject to 


(3-2) 


k(x°)+£ 


dh- 


^ r)x 

k=l OX k 

k= 1 uX k 


Ax k =e t (; i=l,...,p ) 
Ax k <bj (j = 


where X° represents the current design vector and Ar A , is the difference between two 
successive values for the kth design variable. 

In this way, the nonlinear constrained optimization problem is transformed into a linear one. 
There are many methods to solve this linear inequality system, such as simplex method, 
which is used in this report. 

3.1.2.Sensitivity analysis 

For SLP, a major part of analysis cost is sensitivity analysis. In the cureent analysis, the 
structure is analyzed by displacement-based finite element method and thus all the 
sensitivities can be obtained by calculating the difference of displacements before and after 
the change of every shape design variable. Therefore, it is necessary to calculate the 
displacement change with a shape design variable being changed. To avoid repeating the 
time-consuming solution of stiffness equation, the concept of force derivative is employed to 
obtain all the displacement sensitivities by one single analysis, which can be found in [10], as 
follows. 

In a finite element analysis, we have st if fness equation 


Kq = P (3-3) 

where K is the assembled stiffness matrix, q is the vector of displacement parameters and P is 
a (nodal) force vector (vectors). 

For a small change Ax of shape design variable x h we have 
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(K + AK)(q + A q) = P + AP 


(3-4) 


Substituting equation (3-3) into equation (3-4), we obtain 

KAq = A P- A Kq - A kAq (3-5) 


Ignoring the higher order term in equation (3-5) and dividing by Ar on both sides yields 


A q _K \AP-AKq) 

Ax { Ax t 


(3-6) 


The term on the right hand side of the above equation is known as the force derivative of 
shape design variable x t . The differentiation of stiffness and load in equation (3-6) can be 
calculated by 


A K = \B T DBdQ.~ \B T DBdQ. 

JO, JQ. l 

AP = f N T pclPl - f N T pel®. 

JQ. 2 


(3-7) 


where and Q. 2 represent the domains before and after the variation of the design variable 

Xi (see Figure 3-1). B and D are the geometrical matrix and elastic matrix of an element, 
respectively. 



Figure 3-1 Effect of changing a design variable on the boundary of elements 


To get more accurate stress sensitivities, we use the displacement change directly instead of 
the force derivative in this report, which is given by 


da _ a(q + Aq)-a(q) 
dx i Ax t 


(3-8) 


The numerical examples show that the above equation can provide an improved stress 
sensitivity analysis. 

3.1.3 Objective function 

High stress concentration occurs at the boundary of the adhesive layer and it is essential to 
minimize the peak stresses in the adhesive layer by shape optimization. Then, the maximum 
stress, which determines the strength of the adhesive layer, should be selected as the objective 
function. However, the failure modes of the adhesive layer depend on the material property 
and the load applied. Therefore, as a parametric discussion, three kinds of maximum nominal 
stresses are chosen as the objective function in the current analysis, and they are the maximum 
peel stress, the maximum shear stress and the maximum Von-Mises stress, which are defined 
as follows: 


(a) The maximum peel stress a 

(b) The max im um shear stress x = .Jxj + t( 

(c) The maximum Von Mises stress 

^ = ^V (a v -^) 2+(0 . ~ a ,) 2 +( a ; -^v) 2 +6(xJ +T^ t ) 


(3-9) 

(3-10) 


(3-11) 


According to the basic assumption for the adhesive layer, we have cr =cr = x = 0. 
Therefore, the Von Mises stress can be written as 


x = i^2[o;+3(x^+Ti)] 


(3-12) 
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For the sake of convenient, an equivalent stress is used in the analysis 


% = V a v/ 3 + (<+<) (3-13) 

3.1.4 Shape design variables 

Figure 3-2 depicts one quarter of a symmetric bonded patch that is mapped into x-y plane. 
Some points along the periphery of the bonded patch are selected to define the boundary 
shape of the patch. The movements of these points are assumed to be in both the x- and y- 
direction. At the point where v=0 or y=0, the movement is limited in one single direction, i.e. 
j’-dircction or x -direction only. 



Fig.3-2 shape design variables 

At certain stage, the shape change with the movement of boundary point is very small and 
this movement should be neglected in the current iteration. In Figure 3-2, the movement at 
point A in x-direction can only slightly change the shape of the patch so that it can be ignored. 
Therefore, only the movement in y -direction is considered at this point in the current 
optimization step. The magnitude of shape change can be evaluated by calculating the area 
change of the patch in the analysis. 

It is worth noting that all the movements of shape design variables must take place in the 
curved surface. In this analysis, we map the patch into a flat plane and work out the 
magnitude and direction of shape change. Then the patch is mapped back into the curved 
surface. 
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3.1.5. Constraints 


The constraints for shape optimization of bonded patch mainly concern with restraints on the 
objective function, move limits and other design requirement, such as weight and basic shape 
constraints. 


Firstly, the objective function should not be increased for each optimization step. In other 
words, the maximum nominal stress in the adhesive layer must be smaller than that before a 
shape change is introduced. 

Based on the sensitivity analysis described in Section 3.1.2, we have the stress gradients for 
every design variable. Then the stress change at a point with a shape change can be 
approximately evaluated by 

n 

Aa a) = £Vct,; 4) (xf) (i = 1,...,/«) (3-14) 

;=i 


Where Va| ; (1 is the sensitivity vector of the /th design variable at the /th node in the kth 

(current) optimization step, x- ] is the /th design variable vector in the current iteration, and 

x'.' 11 the /th design variable vector in the last optimization step, m is the number of nodes 

where the maximum nominal stress in the adhesive layer occurs, n is the number of the shape 
design variables. 

So the stress constraints can be expressed as 


a a-u +Aa u) =a (*-D + yva w (i ( . 

I l l IJ X J 


Jk- 1 ) 


)<cr 


(k-i) 


(i = l,...,m) 


(3-15) 


where cr ' k 11 denotes the nominal stress on the /th node in the last optimization step. ’ is 
the maximum nominal stress in the last optimization step. 
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According to the current knowledge of stress distribution in the adhesive layer, stress 
concentration only occurs along the periphery of the bonded patch So only the stresses on 
those nodes along the periphery of the bonded patch are considered in equation (3-15). The 
will enable us to employ a much less number of inequality constraints for the problem 
considered. 


The change of the maximum nominal stress with a shape change, i.e., the objective function, 
can be accordingly given by 


Act 


(*) 

max 


" da ik) 

\ ' j max 

J =1 7 


( X J 


(k) 



(3-16) 


where a^ n j ax is the maximum nominal adhesive stress with the jth shape design variable 
x tk " being moved to x'- " + Ax.. Here Ax . is the small change of shape design variable 
used in the sensitivity analysis. 

Secondly, a set of constraints is imposed on the range of Ax admissible in any iteration so 
that every shape change can be kept in a linear range. The choice of the move limits is very 
important to economy and success of the optimization solution. If the move limits are made 
too small, the convergence will be very slow; if they are too large, oscillation may occur. The 
move lim its can be given as 


Xk-l) 


■ d < xj 


<*> 


< x 


<*-n 


+ d 


(3-17) 


where d is the move lim its. 


From the viewpoint of engineering design, the weight of the structure is an important factor 
for structural design. In general, we need to keep the weight in a certain range. For current 
shape optimization, the weight control is implemented by applying a limit on the area of the 
bonded patch. On the other hand, the shape of the patch may reduce infinitely in some 
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direction to become a very narrow strip. To avoid this bad shape, it is necessary to define a 
basic domain that the shape change must take place out of this domain (see Figure 3-2). Then 
we also have some extra structural constraints as follows 


A, < A < A 

l u 


Xj > X 


bj 


(3-18) 


where A is the current area of the patch, A„ and A; are the upper and lower limit of area of the 
patch, respectively, xbj is the minimum length of design variable xj according to the limit of 
basic domain. 


3.2. Adaptive mesh generation 

It is inevitable that mesh topology varies during the optimization process due to a large shape 
change of the bonded patch. To avoid excessive element distortions and unreasonable large 
aspect ratios, even unreal element, an adaptive mesh generation algorithm should be 
employed to generate quadrilateral element. 

In the current analysis, we assume that the curved surface is a cylindrical surface. Then a 
two-dimensional mesh generator is suitable for the requirements. 

The basic idea of adaptive mesh generation is to map the whole structure, including the 
cylindrical surface, the bonded patch, and the adhesive layer, into a flat plane. Then an 
automatic, two-dimensional finite element mesh generator using quadrilateral elements is 
employed to obtain the mesh topology. Finally, the mesh in the flat plane is projected back to 
the cylindrical surface. Because the cylindrical surface can be easily mapped into a flat plane, 
it is shown that the present mesh generation scheme is very efficient and effective. 

3.2.1. Automatic quadrilateral mesh generation [11] 

The finite element method requires that the domain of interest be subdivided into a mesh of 
discrete elements. The accuracy and the expense of the calculations are strongly affected by 
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the ‘goodness’ of the underlying mesh. An ideal mesh is characterized by small elements 
where stress gradients are high, larger elements where stress gradients are small and elements 
that are regular or undistorted in shape. 

The mesher used here is based on the ‘looping algorithm’ first described by Schoofs et al [11] 
and improved by Talbert and Parkinson [12]. The detailed algorithm for this method can be 
found in [12]. We just give a brief description here. 

In this algorithm, the boundary geometry is defined in a continuous closed loop at first. Then 
the size of element on the boundary is specified as per preferred refinement at different 
locations throughout the mesh. According to a definition of visible nodes (to avoid the 
splitting line falling outside the concave loops), the best splitting line linking two visible 
nodes in a closed loop is determined by (see Figure 3-3) 

SL = C|0 + c 2 T + c 3 E (3-19) 

where 0 is the deviation of the split angles from the ideal (tc/2 ), L is the length of the 
splitting line and 8 is the error in element size natching (node placement by predefined 
element size gradient and refinement coefficient). c 1 ,c 2 ,c 3 are constants used as weighting 
values and the following values were proved to give the best results [11] 

c 1 = 0.5, c 2 = 0.3, c 3 = 0.2 (3-20) 

Based on equation (3-19) and (3-20), the original boundary is split recursively until all loops 
contain either four or six nodes. A six-node loop closure is processed separately. 

Finally, the Laplacian method is used to smooth the mesh generated, which can significantly 
improve the distorted elements. 
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Figure 3-3 Splitting line 


3.2.2 Creating a layer of offset elements along the periphery of the bonded patch 

It is well known that high stress concentration occurs along the periphery of the bonded patch 
Therefore, to minimize the stress concentration in this region, it is necessary to have a couple 
of layers of similar offset elements along the boundary of the bonded patch, which should be 
as square as possible, to improve the accuracy of stress analysis. Furthermore, the size of 
these offset elements should be kept as close to a constant as possible during optimization. 

The offset element defined in this way will benefit the numerical analysis in three ways. 
Firstly, the distance from the boundary of the adhesive layer to the Gaussian point, where we 
can evaluate the adhesive stresses, only changes slightly for different shape of the bonded 
patch. Then it is reasonable to compare the stresses on the Gaussian point before and after 
shape change of the bonded patch. Secondly, it can reduce the nodal stress error caused by 
the change of element size. Thirdly, the total number of nodes will not change dramatically 
because the number of automatically generated nodes depends on the size of element 
predefined on the boundary by user. Thus it avoids generating too many nodes, which may 
result in an unacceptable long computing time for the solution of st if fness equation. 

An approach to generate a layer of offset elements has already been presented in [12], which 
is listed in Table 3-1 as case 1-3. However, this method is based on the element size being 

small enough to avoid some special cases, such as two adjacent interior angles are smaller 
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than 120°. In the current analysis, the boundary of bonded patch changes arbitrarily with the 
shape evolution and it is inevitable to have these cases occur in order to keep a nearly 
constant element size along the boundary of patch. Therefore, some extra mesh generation 
methods are presented in this analysis to fulfill the requirement of the stress analysis, which 
are given in Table 3-1 as case 4-6. Figure 3-4 depicts the detailed configuration for all 6 cases. 

The method described above can successfully generate a layer of quadrilateral offset elements 
along the boundary of the patch. But it may fail to generate one more layer of offset elements 
because the element sizes are changed on the new boundary. However, we can interpolate the 
generated offset elements to two layers of offset elements with aspect ratio 1:2, as shown in 
Figure 3-5. 


Table 3-1 Offset node generation schemes 


Case 

Interior angle 

How the offset node is generated 

1 

120° <oc<240° 

Along a bisecting line of the interior angle 

at a distance equal to the element size at 

that node 

2 

a <120° 

As the resultant of the two vectors which 

point to the two adjoining nodes 

3 

a> 240° 

Same as case 2 except the negative of the 

resultant is used 

4 

a, < 120°, a 2 < 120° 

a, ,a 2 are two adjacent interior angles 

Compare the two angles, the small one is 

processed as same as case 2, while the big 

one is processed as case 1 

5 

a, > 240°, a, > 240° 

a, ,oc 2 are two adjacent interior angles 

All two angles are processed as same as 

case 1 

6 

a, <120°, 120° <a 2 <240°,a, <120° 

a, .a 2 ,a are three continuous angles 

Along a bisecting line of a 2 at a distance 

equal to the element size at that node 


16 

























(b) a <120° 


(e) a, > 240°,a, > 240° 



(c) a>120° (f) a, < 120°,a 2 < 120° 

Figure 3-4 Offset node generation method 



Offset elements Interplated offset elements 

Figure 3-5 Interpolating the offset elements 
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3.2.3 Improvement of the distorted elements 


Theoretically, the automatic mesh generation algorithm mentioned above can avoid the 
distorted elements by mesh smoothing if enough refined element sizes are prescribed on the 
boundary of structure. However, it also means that a large number of nodes are generated and 
follows a much longer calculating time. To overcome this drawback, we present a mesh 
improvement scheme. 


As shown in Figure 3-6, there is an element with one of its interior angle (>150°) being much 
larger than other three interior angles after mesh smoothing, which is unacceptable for finite 
element analysis. We can add one more node on its longer side (not including two sides 
forming the larger interior angle). Next, this element and its adjacent elements are refined into 
three elements. Then the mesh is smoothed again. Repeating this procedure for 2-3 times can 
eliminate the distorted elements. Figure 3-7 gives an example for this improvement with one 
time of modification. 



Figure 3-6 Improvement of the distorted elements 



Figure 3-7 Mesh improvement 
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3.2.4 Mesh scheme 


Based on the two-dimensional automatic mesh generation method described above, we can 
subdivide the cylindrical surface and the bonded patch by mapping them into a flat plane. 
First, two layers of offset elements are generated in the bonded patch along its boundary. Two 
more layers of offset elements are also subdivided outside adjacent to the boundary of the 
patch. Then the structure is divided into two continuous close loops, one is in the region of 
patch, and another is in the remained region. For each loop, the automatic mesh generator is 
applied to obtain quadrilateral elements. Then mesh is mapped back to the cylindrical 
coordinate system. To avoid long and narrow shape of the continuous loop, which makes the 
mesher perform badly [11], the loop of the shell is divided into two loops, as shown in Figure 
3-8. A typical mesh generated by the presented scheme is shown in Figure 3-9. 

It should be pointed out that the resulting bandwidth for a final mesh is very poor. So the 
nodes in the region of the bonded patch and remained region are renumbered according to its 
distance to the central point of the structure, respectively. 



Figure 3-8 Mesh zoning scheme 
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(a) Mesh for circular patch 



(b) Mesh for an optimum patch 


Figure 3-9 Typical mesh used in current analysis 
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3.3 Shape optimization scheme 


Based on the theoretical description in the above three sections, the shape optimization 
scheme is proposed as follows 

(1) Initialization. Input the structural data, material properties, mesh generation 
information (such as element size on the boundary, element size gradient, refinement 
coefficient, etc.), move lim its. 

(2) Generate mesh and conduct stress analysis to obtain the stresses on the boundary of 
the bonded patch. 

(3) Sensitivity analysis. Move all the shape design variables from x ,■ to x, + Ax to obtain 
all the linearised coefficients used in SLP. 

(4) Check the current move limits. If the current limits are small enough, the new 
max im um stress is equal to that in the last iteration step and go to step (8). 

(5) Conduct linear programming analysis to find a feasible solution based on the current 
move limits. If no feasible solution is found, the new maximum stress is equal to that 
in the last iteration step and go to step (8). 

(6) Move the boundary of the boned patch to the new shape. 

(7) Generate mesh using the new shape of the patch and conduct stress analysis to obtain 
the stresses on the boundary of the bonded patch. 

(8) Check the maximum stress. If the new maximum stress is larger than that in the last 
iteration step, reduce the move limits and go to step (4). If the new maximum stress is 
smaller than that in the last iteration step, go to step (3). If the new maximum stress is 
very close to that in the last iteration step, the shape optimization is converged and the 
program is stopped. 

A software called B PATCH-Optimizer has been developed and verified through a number of 
examples. 
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3.4 Numerical examples and discussion 


As shown in Figure 3-10, two circular patches with radius r= 15mm are externally bonded to a 
cylindrical shell. These two patches are located at the middle of the cylinder and only 1/8 of 
structure is considered due to its structural symmetry. The overall length of the shell is 
W=300mm. The external radius of the cylindrical shell is R= 150mm. The thickness of both 
the cylindrical shell and the patches are all 2.0mm. The thickness of adhesive layer is 
1=0. 15mm. An internal pressure p=l .OMPa is applied to the cylindrical shell. 



Figure 3-10 A cylindrical shell bonded tow cirular patches on its external surface 

The shell and the patch are assumed to be metallic with Young’s modulus of £’=70GPa and 
Poisson’s ratio of v=0.3. The adhesive used has a Young’s modulus of £ c =2.4GPa and a 
Poisson’s ratio of v c =0.33. 

Both ends of the cylindrical shells are supported by rigid diaphragms, which allow 
displacements only in v-di recti on of the shell. The area of the bonded patch is limited 
between 3.0r 2 and 6.0r 2 . A 1.2rxl.2r rectangle domain in the center of patch is defined as the 
basic domain. The move lim it d=0.35mm is used according to numerical test. 

3.4.1. Shape evolution in the optimization process 

Figure 3-11 depicts the shape evolution in the whole optimization process. The objective 
function is the maximum Von Mises stress evaluated on the node. To be clear, only a few 
selected optimization steps are plotted here. The final optimum shape step is 33. 
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It can be seen that the length of the bonded patch in y -direction keeps decreasing for all the 
shape evolutions. This is because the peak values of peel and shear adhesive stresses increase 
with the circumferential length of the patch, which can be verified by analyzing the 
rectangular patch with different circumferential length. The results show that the 
circumferential length of the bonded patch is the most important shape parameter to the 
adhesive stress for the investigated structure subjected to an internal pressure. From Figure 3- 
11, we can see that the circumferential length of the optimum patch is close to 0.6r for one 
quarter of the patch. This means that the total circumferential length of the optimum patch is 
1.2r and it is coincide with the size of the prescribed basic domain. Therefore, it is 
predictable that the circumferential length of the bonded patch may become shorter if a 
smaller basic domain is given. 

On the other hand, the length of the bonded patch in ^-direction becomes shorter at first and 
then longer. According to the analysis of rectangular patch with different axial length, the 
axial length can slightly affect the peak value of the adhesive peel and shear stresses. 
Obviously, the effect of the axial length on the adhesive stresses changes with the 
circumferential length. 

Table 3-2 gives the maximum adhesive stresses for different shape of patch. Comparing the 
first three rows, it can be found that the comer cut-off cannot decrease the peak adhesive 
stresses. In contrast, the peak adhesive stresses increase very much with the larger corner cut¬ 
off. Furthermore, the optimum shape of the bonded patch in Figure 3-11 shows that the comer 
should be moved outward instead of inward. 

Comparing the patch shapes in step 5 and 10 in Figure 3-11, the presented optimization 
method can adaptively move the shape of the patch from a bad one (step 5) to a good one 
(step 10). In the final several steps, the move of shape design variables is smaller than that in 
other steps. Then the shape change is much smaller than that in other steps, which can be seen 
from step 25 to the final step (8 steps). This shows that the present optimization method can 
precisely evolve into the optimum shape. Therefore, the present method is suitable for shape 
optimization of a patch bonded to a curved surface. 
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3.4.2. Stresses on the node and Gaussian point 


Two optimum shapes of the bonded patch, based on the stress on the node and Gaussian point 
respectively, are shown in Figure 3-12. The objective function is the maximum Von Mises 
stress. The optimum maximum adhesive stresses are given in Table 3-3. It can be found that 
stresses on the node and Gaussian point result in the similar optimum shapes. The Gaussian 
point is just close to the boundary of the patch. So the stresses on the Gaussian point cannot 
provide us enough accurate degree of stress concentration on the boundary of the adhesive 
layer. Consequently, the optimum stresses obtained from the Gaussian stress are slightly 
higher than those from the nodal stress (see Table 3-3). Therefore, all the other results 
discussed in the following are obtained according to the nodal stress. 

3.4.3. Different objective functions 

In Section 3.1.3, we presented three different objective functions, which are the maximum 
peel stress, the maximum shear stress and the maximum Von Mises stress. The purpose to 
define the objective function in this way is to allow consideration of different possible failure 
modes of adhesive layer. 

Figure 3-13 depicts the optimum shapes of the bonded patch by minimizing different 
maximum adhesive stresses. The peak adhesive stresses obtained from these optimum shapes 
are given in Table 3-2. It should be pointed out that the optimum shape of the bonded patch 
obtained from the maximum Von Mises stress as the objective function is used as the initial 
shape of the bonded when we conduct other two analysis. It is clearly shown in Figure 3-13 
that the optimum shapes for different objective functions are similar to each other. Therefore, 
it may not be necessary to conduct several optimization analyses to find an optimum shape of 
the bonded patch to cover all the possible failure modes. Certainly, we must determine which 
stress is the key parameter to the strength of the adhesive layer. 

From Table 3-2, we can see that the peak stresses in the adhesive layer can be reduced 
dramatically by optimizing the geometrical shape of the bonded patch. Especially for the peel 
stress, it can be reduced in a large amount up to 58.8% in this example. In fact, the peel stress 
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is closely related to the circumferential length of the bonded patch. For the limitation of the 
basic domain, the circumferential length of the patch attains the value as per the basic domain. 
The peel stress should be reduced further if a smaller basic domain is given. The shear stress 
and Von Mises stress are reduced by 24.3% and 29.0%, respectively. This shows that the 
shape optimization can efficiently reduce the peak value of adhesive stresses. The objective 
function should be chosen according to the failure mode in its application. 

3.4.4. Effect of curvature 

Figure 3-14 gives the optimum shapes of the bonded patch for different radius of the 
cylindrical shell with the maximum Von Mises adhesive stress as the objective function. The 
optimum adhesive stresses are given in Table 3-4. It can be seen that all the optimum shapes 
have almost same shape in the circumferential direction. The length of the patch in the 
longitudinal direction increases with the radius of the cylindrical shell. This is because the 
length of the patch in the circumferential direction controls the maximum adhesive stress in 
this circumstance. Using as shorter circular length of the patch as possible can improve the 
strength of the adhesive layer. 

3.4.5 Composite patch 

A composite patch bonded to a cylindrical surface is investigated by the present optimization 
method. The composite material is assumed to have following properties 

£j = m.OGPa E 2 = E 3 = \0.3GPa G l2 = G 13 = l.YIGPa 

G 23 = 4.02GPa jj . I2 = jj , 13 = 0.28 ja 23 = 0.28 

Two kinds of stacking sequences are analyzed, which are [0 2 /90 2 /0 2 /90 2 ] s , 
[90 2 /0 2 /90 2 /0 2 ] s . The thickness of each single ply is 0.127mm. Because of the cross-ply 
stacking, we can use the same mesh scheme as metal patch in the analysis. 

Figure 3-15 depicts the optimum shape of the composite patch. The objective function is the 
maximum Von Mises adhesive stress. The stresses for circular patch and optimum patch are 
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given in Table 3-5. It can be seen that there isn’t much difference of the optimum shape for 
two kinds composite patch. Furthermore, they are similar to the optimum shape obtained 
from the metal patch. This is because the composite patch used here is a quasi-isotropic 
laminate. 

Although similar optimum shapes are obtained for these two composite patches, the peel 
stresses in the adhesive layer are very different. As shown in Table 3-5, the peel stress of 
[90 2 /0 2 /90 2 /0 2 ] s patch is much smaller than that of [0 2 /90 2 /0 2 /90 2 ] s patch for both 
optimum and circular patches. This is believed to be mainly caused by the mismatch of 
stiffness in the circumferential direction between the adhesive layer and lower surface ply of 
the bonded patch. For the load case considered here, the patch is subjected to a high stress in 
the circumferential direction [6]. For [90 2 /0 2 /90 2 /0 2 ] s , the first ply adjacent to the 
adhesive layer is 90°, which is in the circumferential direction. The 90° surface ply is less 
stiffer than the aluminum parent shell in the circumferential direction, thus results in a low 
peel stress. In contrast, the peel stress for [90 2 / 0 2 / 90 2 / 0 2 ] s patch is higher than that for 
[0 2 / 90 2 / 0 2 / 90 2 ] 5 patch due to the fact that the CP surface ply is much stiffer that the 
aluminum parent shell. 

The results in Table 3-5 also show that the optimum shape of patch can decrease the stress in 
the adhesive layer significantly. 

3.4.6. Optimization based on thermal stress 

To reduce the thermal residual stress cased by temperature change during the bonding 
process, shape optimization based on thermal stress is performed here. All the geometrical 
sizes of the structure are as same as those described at the beginning of this section, together 
with the same material of cylindrical shell and adhesive layer. The patch is a composite one 
with the same material prosperities as those in Section 3.4.5. The stacking sequence is 
[90 2 /0 2 /90 2 /0 2 ] s . 
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The coefficient of thermal expansion of the metallic shell is assumed to be a = 0.24x10 4 . 
The coefficients of thermal expansion of the composite patch are taken as 

a, = 0.63xl(T 7 , a, = 0.288 x 10 4 , in which subscript 1 indicates the fiber direction and 
subscript 2 the transverse direction. It should be pointed out that the thermal effect of the 
adhesive layer is not considered here. 

The applied thermal load is a temperature change from 220° C down to 20° C. The boundary 
condition is the same as above. 

Figure 3-16 depicts the similar optimum shape of the composite patch. However, the adhesive 
stress remains very high. Table 3-6 gives the adhesive stresses before and after shape 
optimization. The maximum Von Mises stress is reduced by 24.2%. However, the maximum 
peel stress is increased by 5.7%. 

3.4.7. Different boundary support 

All the examples above have same boundary support. To investigate the boundary effect on 
the optimum shape of patch, the both ends of cylindrical shell are fixed. Two kinds of patches 
are considered, which are metallic and composite patches. The stacking sequence of 
composite patch is [902/02/902/02]s- All other parameters are the same as described above. 

Figure 3-17 depicts the shape evolution of metallic patch boned to a cylindrical shell with the 
maximum Von Mises stress as the objective function. A total of 19 steps of shape change 
were required to obtain the optimum shape. Obviously, the circumferential length of patch for 
the optimum shape is much longer than that in Figure 3-11. The width of the patch for the 
optimum shape in the longitudinal direction remains not significantly changed from the radius 
of the initial shape. This shows that the axial constraint can dramatically affect the optimum 
shape. It is believed that the internal pressure and the fixed end constraints create membrane 
stress in the longitudinal direction, which in turn yields high peel stress along the longitudinal 
boundary of the patch. The created high peel stress restraints outward boundary movement in 
the longitudinal direction. 
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Figure 3-18 shows two optimum shapes of metallic patch with the objective function being 
the maximum peel stress and the maximum Von Mises stress, respectively. The two optimum 
shapes almost coincide to each other. 

Figure 3-19 and 3-20 depict the shape evolution and optimum shapes for the composite patch. 
Comparing the optimum shapes with those in Figure 3-15, the difference is similar to that 
observed for the cases of metallic patch. However, it should be pointed out that the optimum 
patches evolved based on different objective functions, i.e., the maximum peel stress and the 
maximum Von Mises stress, are quite different, although their shapes are similar. 

Table 3-7 gives the optimum maximum stress in adhesive layer for metallic and composite 
patch bonded to cylindrical shell with fixed boundary support. Comparing these results with 
those obtained previously, we can find that the maximum Von Mises stress changes slightly. 
However, for the composite patch with the maximum peel stress being the objective function, 
the optimum maximum peel stress decreases significantly, up to 57.0%. For the case with the 
maximum Von Mises stress being the objective function, the maximum peel stress remains 
very high. Thus it is appropriate to choose the maximum peel stress as the objective function 
in this circumstance if the maximum peel stress causes adhesive failure. 


Table 3-2 The maximum adhesive stresses with different shapes of the boned patch 


(The radius of the cylindrical shell = 150mm) 


Shape of the 
bonded patch 

Peel stress 
(MPa) 

Shear stress 
(MPa) 

Von Mises 
stress (MPa) 

K ** 

Square 

14.22 

(14.00,15.00)* 

18.12 

(14.00,15.00) 

19.90 

(14.00,15.00) 

Square with comer 
cut-off 

13.89 

(8.57, 14.92) 

18.96 

(9.61, 14.70) 

20.45 

(8.57, 14.92) 

Circle 

16.38 

(0.00,15.00) 

19.66 

(0.00,15.00) 

21.82 

(0.00,15.00) 

Optimum 

Objective: 

Peel 

6.75 

(19.90, 8.42) 

15.40 

(19.90, 8.42) 

15.74 

(19.90, 8.42) 

58.8% 

Objective: 

Shear 

6.90 

(18.90, 8.40) 

14.88 

(18.90, 8.40) 

15.58 

(18.90, 8.40) 

24.3% 

Objective: 
Von Mises 

7.81 

(13.98, 9.09) 

14.99 

(13.98, 9.09) 

15.49 

(13.98, 9.09) 

29.0% 


*The values in the brace are the nodal coordinate where the maximum stress is located. 


<7 ( is the value of objective function and ct 0 is the corresponding value of the circular patch. 
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Table 3-3 The optimum max im um adhesive stress based on nodal and Gaussian stresses 
(The radius of the cylindrical shell = 150mm) 


Objective function 

Nodal stress (MPa) 

Gaussian Stress (MPa) 

Peel stress 

7.81 (13.98,9.09)* 

8.42 (0.00, 10.32) 

Shear stress 

14.99 (13.98, 9.09) 

15.57 (0.00, 10.32) 

Von Mises Stress 

15.49 (13.98, 9.09) 

16.31 (0.00, 10.32) 

* The values in the brace are the nodal coordinate where tl 

le maximum stress is located. 


Table 3-4 The optimum maximum adhesive stresses for different radius of cylindrical shell 


Radius of Shell (mm) 

Peel stress (MPa) 

Shear stress (MPa) 

Von Mises stress (MPa) 

150 

16.38 (7.81)* 

19.66(14.99) 

21.82(15.49) 

200 

21.45 (9.84) 

26.11 (20.24) 

28.90 (20.91) 

250 

26.59 (12.27) 

32.60 (25.23) 

36.03 (26.70) 

300 

31.61 (15.17) 

39.02 (30.69) 

43.08 (31.35) 


* The values in the brace are the corresponding maximum stress for the optimum patch with 
initial shape of patch being circle. 


Table 3-5 The optimum maximum adhesive stresses for composite patch (Objective function: 


Maximum Von Mises stress on node. The radius of the cylindrical shell = 150mm) 


Stacking sequence 

Patch 

shape 

Maximum 

Peel stress 

Maximum 

Shear stress 

Maximum 

Von Mises 

stress 

(<T, ** 

[90 2 /0 2 /90 2 /0 2 ] 5 

Circle 

12.91 

(0.00,15.00)* 

21.83 

(0.00,15.00) 

23.07 

(0.00,15.00) 

24.1% 

Optimum 

3.72 

(0.0, 9.39) 

17.40 

(13.98,9.00) 

17.51 

(13.98,9.00) 

[0 2 /90 2 /0 2 /90 2 ] 5 

Circle 

18.52 

(0.00,15.00) 

20.42 

(0.00,15.00) 

23.05 

(0.00,15.00) 

28.2% 

Optimum 

8.79 

(2.99,9.31) 

15.78 

(13.97.9.00) 

16.54 

(13.97.9.00) 


*The values in the brace are the nodal coordinate where the max im um stress is located. 


<j a is the value of objective function and <J 0 is the corresponding value of the circular patch. 
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Table 3-6 The optimum maximum adhesive stresses for composite patch (Objective function: 


Maximum Von Mises stress on node. The radius of the cylindrical shell = 150mm) 


Patch 

shape 

Maximum Peel 

stress (MPa) 

Maximum Shear 

stress (MPa) 

Maximum Von Mises 

stress (MPa) 

( CT o -cuM. ** 

Circle 

135.27 

(14.97, 0.98)* 

231.58 

(0.00, 15.00) 

244.42 

(0.00, 15.00) 

24.2% 

Optimum 

142.99 

(21.36,5.89) 

184.09 

(14.00, 9.00) 

185.37 

(14.00, 9.00) 


*The values in the brace are the nodal coordinate where the maximum stress is located. 


<7 ( is the value of objective function and cr 0 is the corresponding value of the circular patch. 


Table 3-7 The optimum maximum adhesive stresses for the fixed end of cylindrical shell 
bonded metallic or composite patch (The radius of the cylindrical shell = 150mm) 


Lay¬ 

up 

Shape of the 

bonded patch 

Peel stress 

(MPa) 

Shear stress 

(MPa) 

Von Mises 

stress (MPa) 

(<T) -COM, ** 

•g 

Circle 

13.06 

(0.00, 15.00) 

17.93 

(0.98, 14.97) 

19.42 

(0.00, 15.00) 


03 

Oh 

O 

OJ 

Optimum 

Objective: 

Peel 

7.75 

(0.00, 12.66) 

15.35 

(0.00, 12.66) 

16.02 

(0.00, 12.66) 

40.7% 

s 

Objective: 

Von Mises 

7.89 

(0.00, 12.66) 

15.32 

(0.00, 12.66) 

15.99 

(0.00, 12.66) 

17.7% 

43 

8 

Circle 

9.94 

(0.00, 15.00) 

19.28 

(0.00, 15.00) 

20.11 

(0.00, 15.00) 


B 

to 

O 

& 

i 

Objective: 

Peel 

4.27 

(9.85, 11.14) 

16.89 

(8.86, 11.29) 

17.09 

(8.86, 11.29) 

57.0% 

5-h 

o 

u 

t 

Objective: 

Von Mises 

7.09 

(12.89, 1124) 

16.82 

(0.00, 12.75) 

17.07 

(0.00, 12.75) 

15.1% 


*The va] 


ues in the brace are the nodal coordinate where the maximum stress is located. 


<J a is the value of objective function and cr 0 is the corresponding value of the circular patch. 
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Figure 3-11 Shape evaluation in the optimization process 
(The objective function is nodal Von-Mises stress) 



Figure 3-12 Optimum shape based on Gaussian and nodal stress 
(The objective function is Von-Mises stress) 
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x (mm) 

Figure 3-13 Optimum shape of patch based on different objective functions 



x (mm) 

Figure 3-14 Optimum shape of patch bonded to shell with different curvature 
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Figure 3-15 Optimum shape of composite patch with different stacking sequence 



Figure 3-16 Optimum shape of composite patch subjected to temperature change 
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Figure 3-17 Shape evolution for the fixed end of cylindrical shell boned metallic patch 



Figure 3-18 Optimum shape based on different objective functions for the fixed end of 

cylindrical shell boned metallic patch 


34 



16 



Figure 3-19 Shape evolution for the fixed end of cylindrical shell boned composite patch 



Figure 3-20 Optimum shape based on different objective functions for the fixed end of 
cylindrical shell boned composite patch 
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4. Optimization of the ply drop-offs of composite patch 


The stress concentration, which occurs in the adhesive layer along the periphery of the bonded 
patch, can dramatically reduce the bonding strength. It is thus desirable to minimize the 
adhesive stresses by optimizing the shape of the curved bonded patch. However, even for the 
optimum shape of the patch, the stress concentration near the boundary of the adhesive layer, 
which is caused by the discontinuities of material and geometry on the bonding interface, can 
still be high too. Therefore, it is necessary to relieve the stress concentration on the boundary 
of the adhesive layer by employing other techniques. 

There are two methods available to improve the stress concentration on the interface between 
the different materials. From the analysis of the bi-material connection, it was found that the 
stress singularity might disappear by choosing appropriate combinations of material constants 
in the two materials when the geometry is fixed [13, 14]. In the second method, changing the 
geometric shapes, such as tapering, along the periphery of the bonding interface can also 
reduce the stress concentration, according to the research results for flat adhesive joints [6, 15]. 

In general, the material and the geometry of the curved surface are pre-determined before we 
start to optimally design the bonded patch. Then only the material and geometry of the bonded 
patch can be select as design variables. This makes it extremely difficult to choose an 
adequate combination of material constants for the bonded patch because the materials 
available are in general limited. Composite material provides us more opportunities to design 
the properties of material by designing its lay-up. This method will be discussed in the next 
section. In this section, we focus on optimum design of the geometry of the bonding interface 
in the through thickness direction, i.e., optimum tapering. 

Based on the current knowledge of steel wedging and flat bonding design, the wedge angle in 
the both adherents is a key geometric parameter for the alleviation or ever disappearance of 
stress concentration. The stress concentration can even be removed by cutting the adherent at 
an appropriate wedge angle on its boundary. Moreover, a wedge angle is very easy to 
implement for composite patch by the ply drop-off. Therefore, the wedge or tapering angle is 
selected as the optimization parameter. 
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The wedge angle is suitable for flat joint or repair design because its cutting plane is a flat one. 
This cutting plane becomes a curved one if the wedge is applied to a curved joint design. For 
the sake of convenience, we use a similar parameter, the drop-off width, to replace the wedge 
angle, which is described in the following sections. 

In this section, a detailed stress analysis method of the bonded composite patch with the ply 
drop-offs is presented first. Then a mesh scheme for the ply drop-off is given to conduct the 
stress analysis. It should be noted that the ply drop-off in the current analysis is limited to be 
uniform to reduce the computational time. As a matter of fact, the present method can also be 
applied to conduct an analysis of the non-uniform ply drop-off. Finally, the effects of ply 
drop-off on the adhesive stresses for the different shape of patch are investigated. 

4.1. Definition of the ply drop-off width 

Figure 41 depicts a curved composite patch with the ply drop-off is bonded to a curved 
surface. The ply drop-off takes place on the all the boundary of every ply of the bonded patch. 
The definition of the wedge angle and the cutting plane are shown in Figure 42. It can be 
clearly seen that the wedge angle can fully describe the ply drop-offs for a flat bonded patch 
(see Figure 4-2(a)). Accordingly, when the geometry of flat is projected to a curved plane, (as 
shown in Figure 4-2(b)), the flat cutting plane becomes a curved one. Then the wedge angle is 
unsuitable to define the ply drop-offs for the curved patch. Therefore, the drop-off width d is 
introduced to define how every ply is dropped off (see Figure 4-2(b)). 



Figure 4-1 A composite patch with the ply drop-offs is boned to a curved surface 
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(a) Flat patch 



Figure 4-2 The wedge angle and cutting plane of patch with the ply drop-offs 

4.1.1 Ply drop-off scheme for special shape of patch 

For some special shapes of patch, such as a circular or square patch, we can define the drop¬ 
off width as the normal distance from the boundary of one ply to that of adjacent ply. It means 
that the drop-off width keeps uniform along all the boundary of the bonded patch. Moreover, 
the drop-off region can be extended to its central point, which enable us design a long enough 
total drop-off width from its external boundary to its central point. Figure 4-3 gives the 
definition of the ply drop-off for these two kinds of typical shape of patch. 




Figure 4-3 Definition of the ply drop-off width 


38 



































4.1.2. Ply drop-off scheme for general shape of patch 


Although it is easy to find the new boundary for a new ply drop-off for the shape of patch 
mentioned above by the definition of ply drop-off width in Section 4.1.1, it may fail to find 
such a new boundary for a new ply drop-off if an arbitrary shape of patch is given. For 
example, for the optimal shapes of the bonded patch, which are obtained from shape 
optimization, only a very small range of drop-off width can be applied to get the new 
boundary with uniform drop-off width. However, a longer drop-off width can reduce the 
highly concentrated adhesive stress on the boundary of the adhesive layer. Therefore, we need 
to design a scheme to calculate the new boundary for each ply drop-off. 

In this report, we only consider the cylindrical surface bonded with a patch. Furthermore, we 
assume that the bonded patch is a symmetrical one, which has two symmetrical axes. Then the 
ply drop-off width is defined as the distance between two boundaries on its symmetrical axes, 
as shown in Figure 4-4. To be convenient, we can calculate the new boundary of the drop-off 
in a flat plane. In Figure 4-4, x and y are the two symmetrical axes of the bonded patch, with 
the length 2a and 2b, respectively. The ply drop-off width cl is defined on these two axes. 
Then the new boundary of the drop-off can be given by 


b-d 

y i =^0 •— 7 — 

b 


(4-1) 


where xo, yo are the coordinate of the node on the initial boundary of patch, jci, yi arc the 
coordinate of the node on the new drop-off boundary. It can be verified that the drop-off 
scheme described here can give the same drop-off shape as that obtained in Section 4.1.1 if 
the initial shape of the bonded patch is circular and square. 

Figure 4-5 depicts an example of drop-off shape of the bonded patch obtained by the method 
presented above. 
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Figure 4-4 Ply drop-off definition for general shape of patch 



Figure 4-5 Ply drop-off of a general shape of patch (Projected into a flat plane) 


4.1.3. Uniform ply drop-off width 

It should be pointed out that it is possible to define different drop-off width between two 
different adjacent plies, which means that the cutting plane at the end of the patch is a curved 
one instead of a flat one even for a flat patch n Figure 42(a). However, to simplify the 
problem and reduce the computing time, a uniform drop-off width is assumed to the bonded 
composite patch. Then we can reduce the design parameter to a single uniform drop-off width, 
which can be found by employing an one-dimensional optimization method. 
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4.2. Structural analysis of the bonded composite patch with the ply drop-offs 

For the shell element used here, all the nodal displacements in an element are referred to the 
mid-plane of the element and the thickness remains constant in the element. In the former 
analysis, the thickness of both the bonded patch and the curved surface are constant and the 
adjacent elements share the same nodes (defined in their mid-planes) in their joining plane. So 
all the nodal displacements are referred to the mid-plane of the bonded patch or the curved 
surface. 

However, the existence of the ply drop-offs makes the thickness of the patch change with the 
ply drop-offs. This results in different mid-planes for any two adjacent shell elements along 
the ply drop-off direction because of both elements having different thickness. It follows that 
the nodes defined in the adjacent two elements along the ply drop-off direction are located at 
the different position in their joining normal, which results in tow nodes in one normal. These 
two nodes should be reduced to one node according to the assumed displacement field. Then 
we can define a common reference plane, called master plane. The node in this master plane is 
called master node. Accordingly, the nodes defined in the adjacent elements are called slave 
node. The master node is used in the total stiffness matrix, which indicates that the total 
stiffness matrix is referred to the master plane. Before assembling the element stiffness matrix 
(referred to its own mid-plane) into the total stiffness matrix, the element stiffness matrix must 
be transformed from the slave node to the master node. It is worthy nothing that the master 
node can be arbitrary selected in the joining normal. 

In Figure 4-6, we define the node in the mid-plane of the bonded patch without the ply drop- 
offs as the master node, the node in mid-plane of the shell element as slave node. All the 
nodal displacements in the patch are referred to these master nodes. It is then convenient to 
transform all the nodal displacements from the mid-plane of an element to the corresponding 
master node. 

According to the first order shear deformation theory for plate and shell, we have the 
relationship of the displacements between the slave node and the master node 
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Figure 4-6 Master and slave nodes in the ply drop-offs 
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where h is the distance between the corresponding pairs of slave and master nodes. The 
subscript ‘ m ’ denotes the displacements of the master node, and the subscript V means the 
displacements of the slave node. 

Rewrite equation (4-2) in a vector form 

{q.} = m{q m } (4-3) 

where {q m } is the nodal displacement vector referred to the mater node, {q s } is the nodal 
displacement vector referred to the slave node. 

Then all the element stiffness matrices are evaluated on their own mid-plane and all the nodal 
displacements are referred to the slave nodes. The strain energy in the element can be given 
by 
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u =^\ A iBtq s } e ) T [D][B\q s YdA 


(4-4) 


where } is the nodal displacement vector in the element referred to the slave node. [5] is 

the element geometrical matrix. [D] is the elastic matrix referred to the mid-plane of the 
element. A is the area of the element. 

Substituting equation (4-3) into equation (4-4), we have 

u =\l{B\\rtq„y) T [D\B\Ttq m YdA (4-5) 

According to the variation theory, we can obtain the element stiffness matrix evaluated on the 
master node 

[/a=£([sId) r MsMM =[[r]W[DM[d/A 

= M r ([ [fit [Dp]4d = [7f [K, ]■ [r] (4 ‘ 6) 

where is the element stiffness matrix evaluated on the master nodes and \K ]' is the 

element st if fness matrix evaluated on the slave node. 

Therefore, the element stiffness matrix is calculated in terms of the displacements referred to 
the mid-plane of element, i.e., the slave nodal displacements. Then the element stiffness 
matrix is transformed into that under master nodes by equation (4-6). Finally, the element 
st if fness matrix is assembled to the global st if fness matrix. 

4.3. Mesh scheme 

According to the structural analysis method described above, every ply drop-off must be 
divided into different elements. To improve the precision of stress analysis, sever similar 
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elements should be added before and after the ply drop-offs. Then the automatic mesh 
generator can be used to generate a mesh for the remaining part of the structure. Figure 4-7 
depicts a mesh scheme of a square bonded patch with 4-layers of the uniform ply drop-offs. 
There is also one layer of similar elements adjacent to the ply drop-off region. In fact, there 
should be more layers of similar elements outside the ply drop-offs in the analysis. Figure 4-8, 
4-9, 4-10 and 4-11 depict typical meshes generated by this mesh scheme. 

The total number of node increases dramatically with the total number of ply of the composite 
patch because every ply drop-off must be divided into at least one element. To improve the 
computational efficiency, the structural parts that don’t have ply drop-offs and are 
automatically meshed can be analyzed first as substructures, which can be incorporated into 
the ply drop-off zone. 



Figure 4-7 Mesh scheme for a square patch with the ply drop-offs (One quarter of stmcture) 
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A-A 


(a) Mesh in a flat plane 



(b) local details of A-A 


Figure 4-8 Mesh generated in a flat plane for a circle patch with ply drop-offs boned to a 

cylindrical shell 
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(a) Global view of structure 



(b) Local details of A-A 


Figure 4-9 Typical mesh for a circle patch with ply drop-offs bonded to a cylindrical shell 
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(a) Global view of structure 



(b) Local details of A-A 


Figure 4-10 Typical mesh for a square patch with ply drop-offs bonded to a cylindrical shell 
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(a) Global view of structure 



(b) Local details of A-A 

Figure 4-11 Typical mesh for the optimum shape of patch with ply drop-offs bonded to a 

cylindrical shell 
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4.4 Optimization method 


Because only one single drop-off width is employed as the design variable, a one-dimensional 
direct search method is employed to find the optimum drop-off width. In this report, the 0.618 
division method is used. 

As shown in Figure 4-12, x is the optimization parameter and fix) is the objective function. It 
is possible that there may exist several peak values in the whole range of x and the maximum 
and minimum values of objective function can be obtained from these peak values. However, 
it is difficult to know how many peak values in the whole range of x. Therefore, the whole 
range of optimization parameter should be divided into several segments. For each segment, 
we can assume that there are only two peak values, which stand for the maximum and 
minimum value respectively. These two peak values can be searched iteratively by using the 
0.618 division method. Then we can obtain the maximum or minimum value of objective 
function by comparing the peak values in these segments. Accordingly, the optimum ply 
drop-offs width is the corresponding width that minimizes the objective function. 



Figure 4-12 One-dimensional direct search scheme 

4.5 Verification of the proposed FEA module 

As shown in Figure 4-13, a cantilever curved beam bonded with a patch-beam is subjected to 
a tangent load T at its free end. The curvature radius of beam is R=\ 50mm and its sector angle 
is 30°. The circumferential length of the bonded patch-beam is 15.0mm. The thickness of the 
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parent beam and the patch-beam are 2.0mm and 1.0mm respectively. The thickness of 
adhesive layer is 0.15mm. The tangent load is equal to 1.0 N/mm. 



Figure 4-13 Cantilever curved beam with drop-offs subjected tangent load at its free end 

Total 8 drop-offs are considered in the patch-beam. The uniform drop-offs width is 0.25mm. 
The circumferential length of drop-off region is 2.0mm. Both the parent beam and the patch- 
beam are assumed to be aluminum with a Young’s modulus of i?=70GPa and a Poisson’s ratio 
of v=0.3. The adhesive has a Young’s modulus of E c =2.4GPa and a Poisson’s ratio of v c =0.33. 

Commercial finite element analysis software STRAND 7 is used to analyze the problem and 
to serve as a benchmark for the proposed FEA module. Plane strain element is used to model 
both the beam (parent beam and patch-beam) and the adhesive layer. For the patch-beam, 8 
elements are divided across the overall thickness and the same size 0.125mm is used in its 
circumferential direction. The same element size is applied to the parent beam. The 
circumferential size of the adhesive layer is 0.125mm too. Only one element is used across the 
adhesive layer’s thickness. For the shell element in the proposed FEA module, the 
circumferential size is equal to the drop-off width. The width of beam is taken as 10mm for 
shell analysis. 

Figure 4-14 depicts the adhesive stresses distribution near the free end of the overlap region 
along the circumferential direction for both with and without drop-off at the end of patch- 
beam. It can be seen that the predictions given by the present FEA module correlate well with 
the results of STRAND7. 
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An inspection of the peak values of the adhesive stresses in Figure 4-14 shows that the peak 
value of the peel stress decreases with the existence of drop-offs in the patch-beam, and the 
absolute value of compressive stress becomes smaller too. However, the peak value of the 
shear stress remains unchanged. The stress distribution patterns remain the same with the peak 
value near the end of the adhesive layer for both peel and shear stress. 



(a) No ply drop-off 



(b) Ply drop-off width=0.25mm 

Figure 4-14 Adhesive tresses distribution along the circumferential direction 
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4.6. Numerical examples 


The same cylindrical structure analyzed in 3.4 is considered to investigate the effect of the ply 
drop-offs of the bonded patch on the peak adhesive stresses. It should be noted that the patch 
is made of composite materials. The bad and boundary condition are as the same as those 
given in Section 3.4. 

The composite material used here is assumed to have the following properties 

E ] = 181.0GPa E 2 =E 3 =10.3GPa G 12 =G 13 =7.17GPa 

G 23 = 4.02GPa jli , 2 = ja 13 = 0.28 = 0.28 

The stacking sequence of the composite patch are [90s], [Os], [90/0/90/0] s and 
[0/90/0/90] s , respectively. The thickness of a single ply of composite is 0.127mm. Then the 

thickness of the patch is 1.016mm. Because only cross-plies are taken into account, we can 
take advantage of its structural symmetry to save the computational time by analyzing one 
quarter of the structure. 

In the meantime, the uniform ply drop-offs width is defined for each ply drop-off step. To 
investigate the effect of the residual thickness (the thickness of the first drop-off step starting 
from the boundary of patch and adjacent to the adhesive layer), the residual thickness is 

defined as 1, 2 and 3 plies. Then total 7, 6 and 5 steps of ply drop-offs are considered 

accordingly. 

Three shapes of the bonded patch are analyzed here, i.e., circular patch, rectangular patch and 
the optimum shape of patch obtained by the shape optimization analysis in Section 3. For 
circular and rectangular patches, [90s] and [0s] are analyzed. For the optimum shape of patch, 
[90/0/90/0],; and [0/90/0/90 ] s are considered. 

It is obvious that the drop-off width cannot be enlarged unlimitedly because of the limited size 
of the bonded patch. The range of drop-off width is then given to avoid the unrealistic drop- 
offs. In this example, the uniform drop-off width for circular and square patch is from 0.1mm 
to 1.5mm. The total width of drop-off region can be in a range of 0.7mm to 10.5mm. For the 
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optimum shape of patch, the drop-off width is from 0.1mm to 1.0mm with a drop-off region 
ranging from 0.7mm to 7.0mm. 


4.6.1 Circular patch 

Figure 415 to Figure 420 depict the maximum normalized adhesive stresses for circular 
composite patch with different residual thickness with the ply lay-up [0x|. [90s]. The 
maximum adhesive stresses with and without ply drop-off are defined as a max and CT max0 , 
respectively. 

From Figure 415 to 417, it is obvious that the introduction of ply drop-off can reduce the 
stress concentration quite remarkably. The maximum peel stress for [0s] patch with the 
residual thickness 0.127mm is only 21.2% of its initial value when the ply drop-off width is 
1.5mm (see Figure 413). 62.2% and 60.1% of the initial values exist for the maximum shear 
and Von Mises stresses under the same condition. Tie maximum stresses are decreased by 
78.8%, 37.8% and 39.9%, respectively. It is easy to find that the maximum peel stress is 
reduced much more than the maximum shear and Von Mises stress for the same drop-off 
width. In fact, the applied load is mainly transferred from the cylindrical shell to the bonded 
patch by shear stresses in the adhesive layer, and thus the shear stress cannot be reduced 
dramatically. 

However, the maximum stresses don’t continuously decrease with the extension of the ply 
drop-off width. At certain stage, the maximum stresses may increase slightly with an 
increasing drop-off width. For the [90s] patch shown in Figure 418 to 420, the maximum 
peel stress decreases at first (drop-offs width=0.0 - 0.15mm), then it increase with the drop¬ 
off width in the range of 0.15mm to 0.6mm, and eventually the peel stress decrease again. 
This shows that it is necessary to choose proper drop-off width by optimization analysis to 
minimize the peak value of adhesive stresses. 

Comparing the stress level for different residual thickness of the patch, it can be easily found 
that the objective stresses for the thicker residual thickness are generally higher than that for 
thinner residual thickness with the same drop-offs width. This shows that the thinner residual 
thickness results in more reduction in adhesive stress. 
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On the other hand, the effect of ply drop-off on the adhesive stresses changes with the 
stacking sequence of composite patch. For the same drop-off width and same residual 
thickness, the adhesive stresses of [Os] patch are reduced much more than that of [90s] patch. 
As clearly shown in Table 4-1, the adhesive stresses of [Os] patch are much lower than that of 
[90s] patch. It indicates that the effect of ply drop-off for the [Os] patch is more remarkable 
than that for the [90g] patch for this load case. 

It can also be seen that there are some uncontinuous points in these curves in Figure 4-15 to 4- 
20. This is because the position of the maximum adhesive stresses in the bonded patch 
changes with different ply drop-off width. Table 4-2 shows the coordinate of the maximum 
adhesive stresses [Os] patch with different ply drop-off width. 

4.6.2 Square patch 

The results of maximum adhesive stresses for square patch are given from Figure 4-21 to 4-26. 
It can be seen that the maximum peel and Von Mises stress are reduced in similar fashion as 
that of circular patch. For the [Os] patch, the adhesive stresses of square patch are reduced 
slightly less than that of the circular patch with the same drop-offs width and residual 
thickness. However, for the [90s] patch, the adhesive stresses of the square patch are reduced 
slightly more than that of the circular patch with the same drop-off width and residual 
thickness. This shows that the effect of ply drop-off $ affected by the stacking sequence of 
composite patch. 

It should be noted that the maximum shear stress of the [Os] square patch increases with the 
introduction of small ply drop-off (see Figure 4-21 to 4-23). For the residual thickness equal 
to 0.381mm, the maximum shear stress remains larger than that without drop-off. This is 
because we define the maximum shear stress as a combination of two components x , T xz , 
which is given in equation (3-13). The ply drop-off decreases the component x for the 
whole range of drop-off width. But the component x v _ increases rapidly with a very short 
drop-off width, which is caused by the load and boundary condition. Meanwhile, the position 
of the maximum shear stress is different too. The numerical results show that thinner residual 
thickness is preferred. 
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4.6.3 Patch with the optimum shape 


To investigate the effect of ply drop-offs on the patch with the optimum shape obtained by 
shape optimization in Section 3, two types of stacking sequences of composite patches, 

[90 / 0/ 90/ 01 s . and [ 0/ 90/ 0/ 90] v , are analyzed using the presented module. The maximum 

adhesive stresses with different ply drop-off width for these two patches are shown in Figure 
4-27 and 4-28. 


Obviously, the existence of ply drop-off can reduce the maximum adhesive stresses for these 
two cross-ply patches. But it can be seen that the maximum peel stress for [0/90/0/90 ] s 
decreases first, then increases dramatically and finally decreases It is evident that it should be 
careful to introduce drop-off to composite patch bonded to a curved surface. The adhesive 
stresses may increase if the drop-off width is not a adequately selected, which clearly 
demonstrate the necessity of optimal selection of tapering. 


Table 4-1 Maximum adhesive stresses with the optimum drop-off of composite patch 

(residual thickness=0.127mm) 



Lay-up of 
patch 

Drop-off 

width 

Maximum 
peel stress 

Maximum 
shear stress 

Maximum Von 
Mises stress 

Circle 

[0 8 ] 

0.0 

11.37 

(0.98, 14.96) 

9.43 

(9.13, 11.90) 

10.18 

(1.96, 14.87) 

1.5 

2.42 

(0.00, 15.00) 

5.87 

(15.00, 0.00) 

6.11 

(15.00, 0.00) 

[90 8 ] 

0.0 

14.79 

(0.0, 15.00) 

20.06 

(0.50, 14.99) 

21.80 

(0.00, 15.00) 

1.5 

8.23 

(0.00, 15.00) 

17.17 

(0.00, 15.00) 

17.82 

(0.00, 15.00) 

Square 

[0 8 ] 

0.0 

9.47 

(2.00, 15.00) 

7.33 

(13.50, 15.00) 

9.09 

(13.50, 15.00) 

1.5 

2.55 

(14.00, 15.00) 

5.90 

(15.00, 15.00) 

5.91 

(15.00, 15.00) 

[90 8 ] 

0.0 

17.10 

(14.50, 15.00) 

20.61 

(14.00, 15.00) 

22.71 

(14.00, 15.00) 

1.5 

7.34 

(9.00, 15.00) 

15.71 

(9.00, 15.00) 

16.28 

(9.00, 15.00) 

Optimum 

shape 

[90/0/90/0], 

0.0 

9.52 

(16.97, 8.77) 

15.42 

(19.95, 8.48) 

16.01 

(16.97, 8.77) 

1.0 

4.72 

(12.98, 9.02) 

12.94 

(13.98, 9.00) 

13.21 

(13.98, 9.00) 

[0/90/0/90], 

0.0 

12.58 

(13.98,9.00) 

14.61 

(19.95, 8.48) 

15.94 

(16.97, 8.77) 

1.0 

10.49 

(13.98, 9.00) 

10.18 

(13.98. 9.00) 

11.85 

(13.98, 9.00) 
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Table 4-2 Position of max im um adhesive stresses in the [Os] circle patch with residual 

thickness=0.127mm 


Drop- 

offs 

width 

(mm) 

Maximum 

peel stress 
(MPa) 

X 

(mm) 

y 

(mm) 

Maximum 

shear 

stress 

(MPa) 

X 

(mm) 

y 

(mm) 

Maximum 

Von Mises 

stress 

(MPa) 

X 

(mm) 

y 

(mm) 

0.00 

11.37 

0.98 

14.96 

9.43 

9.13 

11.90 

10.18 

1.96 

14.87 

0.10 

8.57 

0.00 

15.00 

7.63 

9.13 

11.90 

8.27 

4.82 

14.20 

0.15 

8.30 

0.00 

15.00 

7.51 

9.13 

11.90 

8.07 

4.82 

14.20 

0.20 

7.75 

0.00 

15.00 

7.38 

9.13 

11.90 

7.82 

4.82 

14.20 

0.25 

7.15 

0.00 

15.00 

7.22 

9.13 

11.90 

7.54 

4.82 

14.20 

0.30 

6.57 

0.00 

15.00 

7.05 

9.13 

11.90 

7.25 

5.74 

13.86 

0.35 

6.07 

0.00 

15.00 

6.88 

9.89 

11.28 

6.99 

5.74 

13.86 

0.40 

5.65 

0.00 

15.00 

6.75 

9.89 

11.28 

6.94 

15.00 

0.00 

0.45 

5.27 

0.00 

15.00 

6.63 

10.61 

10.61 

6.91 

15.00 

0.00 

0.50 

4.91 

0.00 

15.00 

6.51 

11.28 

9.89 

6.88 

15.00 

0.00 

0.55 

4.64 

0.00 

15.00 

6.42 

11.90 

9.13 

6.84 

15.00 

0.00 

0.60 

4.39 

0.00 

15.00 

6.35 

12.47 

8.33 

6.80 

15.00 

0.00 

0.65 

4.16 

0.00 

15.00 

6.32 

15.00 

0.00 

6.78 

15.00 

0.00 

0.70 

3.96 

0.00 

15.00 

6.29 

15.00 

0.00 

6.73 

15.00 

0.00 

0.75 

3.79 

0.00 

15.00 

6.29 

15.00 

0.00 

6.72 

15.00 

0.00 

0.80 

3.63 

0.00 

15.00 

6.26 

15.00 

0.00 

6.68 

15.00 

0.00 

0.85 

3.50 

0.00 

15.00 

6.24 

15.00 

0.00 

6.64 

15.00 

0.00 

0.90 

3.38 

0.00 

15.00 
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Figure 4-15 The maximum stresses for [Os] circular composite patch with the residual 

thickness=0.127mm 



Ply drop-off width (mm) 


Figure 4-16 The max im um stresses for [Os] circular composite patch with the residual 

thickness=0.254mm 
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Figure 4-17 The maximum stresses for [Os] circular composite patch with the residual 

thickness=0.3 81 mm 



Ply drop-off width (mm) 


Figure 4-18 The max im um stresses for [90s] circular composite patch with the residual 

thickness=0.127mm 
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Ply drop-off width (mm) 


Figure 4-19 The max im um stresses for [90s] circular composite patch with the residual 

thickness=0.254mm 



Ply drop-off width (mm) 


Figure 4-20 The max im um stresses for [90s] circular composite patch with the residual 

thickness=0.3 81mm 
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Ply drop-off width (mm) 


Figure 4-21 The maximum stresses for [Os] square composite patch with the residual 

thickness=0.127mm 
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Figure 4-22 The maximum stresses for [Os] square composite patch with the residual 

thickness=0.254mm 
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Figure 4-23 The maximum stresses for [Os] square composite patch with the residual 

thickness=0.381 mm 
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Figure 4-24 The max im um stresses for [90s] square composite patch with the residual 

thickness=0.127mm 


61 




1.2 



C/3 

C/3 


<D 

Vh 

C/3 


0.6 


-d 



o 

23 


0.2 


—*— Maximum peel stress 
—■—Maximum shear stress 
* Maximum Von Mises stress 


0 



0 0.3 0.6 0.9 1.2 1.5 


Ply drop-off width (mm) 


Figure 4-25 The max im um stresses for [90s] square composite patch with the residual 

thickness=0.254mm 
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Figure 4-26 The max im um stresses for [90s] square composite patch with the residual 

thickness=0.381mm 
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Ply drop-off width (mm) 

Figure 4-27 The max im um stresses for optimum shape of [0/90/0/90]s composite patch 
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Figure 4-28 The max im um stresses for optimum shape of [90/0/90/0]s composite patch 
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5. Optimization of the ply lay-ups of a composite bonded patch 

As mentioned in Section 4, the alternative approach to improve the stress concentration on the 
interface between the different materials is to choose an appropriate combination of material 
constants for the two candidate materials when the geometry is fixed. Although this 
combination can be theoretically found, it is not easy to implement because of the limitation 
of available materials. Therefore, this concept is not widely applied to the metal patch. 
However, composite materials enable us to adjust the global proprieties of a laminate 
remarkably by changing its ply lay-ups. So it is practically desirable to search for optimum 
ply lay-ups to reduce the stress concentration in the adhesive layer. 

In the current analysis, the SLP method is employed to conduct the optimization analysis of 
the ply lay-ups of a composite bonded patch. The objective is to minimize the adhesive 
stresses and the design variables are the fiber orientation angles of all plies in the composite 
patch. The shape and thickness of the bonded patch are assumed fixed. The analysis approach 
is presented first and follows some numerical examples. 

5.1. Optimization method 

The SLP method, briefly described in 3.1, is applied to conduct the optimization analysis. 
The formulations used are given in the following. 

5.1.1 Objective function 

The purpose of optimization is to minimize the adhesive stresses. So the adhesive stresses are 
selected as the objective function. Here, the maximum peel stress and Von Mises stress (see 
equation (3-9) and (3-11)) in the adhesive layer are used as the objective function, 
respectively. 

5.1.2 Design variables 


Only the fiber orientation angles of every ply in the composite patch are chosen to be design 
variables and they are assumed to take values of full degree(s) between -90 and 90 degree(s). 
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The shape and thickness of the bonded patch are not design variables and they are assumed to 
be fixed in the current analysis. 

5.1.3 Sensitivity analysis 

The sensitivity analysis for the ply lay-ups optimization is basically similar to that for the 
shape optimization in equation (3-6). The difference is how to calculate the increment of the 
global stiffness and applied load. 

Because only the fiber orientation angles of all plies in the composite patch are the design 
variables, we only need to calculate the increment of the global stiffness and load caused by 
the change of fiber orientation. 

As assumed above, the shape and thickness of the bonded patch are fixed in the whole 
analysis, i.e., the geometry of structure keeps unchanged. This results in that there is no 
change of the applied load. Therefore, we have 

AP = P'—P = 0 (5-1) 

For the change of fiber orientation of each ply in the composite patch, the elastic matrix of 
bonded patch is changed accordingly 

D'= D + AD (5-2) 

where D' is the elastic matrix of bonded patch with new ply-ups. D is the elastic matrix of 
composite patch before the ply orientation angle is changed. AD is the increment elastic 
matrix. 


For the differentiation of stiffness, we have 
A K = \b t D'BcIQ- \b t DBcIQ 

Jq Jq 

= f B T (D + AD)BclQ - f B T DBclQ 

Jq 


(5-3) 
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where A k is the increment global stiffness matrix caused by the change of fiber orientation 
angle. The rest terms are the same as are given in equation (3-7). 

Then the displacements, before and after the fiber orientation is changed, can be obtained as 
follows 


q = K~ l P 

q'=(K') l P'=(K+AK) l P 


(5-4) 


Ignoring the second small term and following the similar procedure from equation (3-4) to (3- 
6), the variation of displacements caused by the change of the fiber orientation angle can be 
given by 

A q = q'-q = -K l AKq (5-5) 

Then substituting equation (5-5) into equation (3-8), we can obtain all the stress sensitivities 
at a node or a Gaussian point for every change of fiber orientation. It is worth pointing out 
that x t stands for the fiber orientation angle of the ith ply in the composite bonded patch here. 


5.1.4 Constraints 

It is obvious that the stress constraints for the optimization of the ply lay-ups are the same as 
those for the shape optimization described in Section 3. Therefore, we have 


^ da 


(*> 


(/ = !,...,/«, j=l, ...,n) 


L -‘ dr 

J =i >X j 


(5-7) 


where n is the number of plies in the bonded patch and x denotes the fiber orientation angle. 
The rest terms are as same as those in equation (3-15). 
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Similarly, the objective function is given by 


(5-8) 

OXj 

For each ply of the composite bonded patch, the fiber orientation angle is limited in a range 
of -90 to +90 degrees because of the periodical characteristics of composite material. Then 
we have 


E 


-90 < x ( p < 90 (j=\,...,n) 


(5-9) 


Finally, a move lim it is imposed to design variables 

-d< x, 1 - xf-" < d (5-10) 

It is worth noting that the optimization scheme presented above may result in an arbitrary 
stacking sequence of composite patch. Then the symmetry across the thickness of optimal 
patch cannot be maintained. However, symmetrical lay-ups of composite patch are generally 
applied in engineering structure. One more constraint, therefore, is imposed for simplicity to 
the optimization method. The corresponding two symmetrical plies (refer to the mid-plane of 
the composite patch) are assumed to change in a same angle in each iteration step. 

5.2 Numerical example 

The cylindrical structure described in Section 3.4 is taken to study optimization of the ply 
lay-ups of composite patch. Both ends of the cylindrical shells are fixed to three translation 
displacements, which allow rotational displacements only. An internal pressure p=\ .OMPa is 
applied to the cylindrical shell. 

The same material is applied to the cylindrical shell and the adhesive layer with thickness 
2.0mm and 0.15mm respectively. Other size of cylindrical shell is as same as that given in 3.4. 
The composite patch has the same material properties as given in Section 4.6. 
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Two kinds of shape of patch are considered here, which are circle and square patches. The 
initial stacking sequences of the bonded patch are [90i6] and [90s]. The thickness of a single 
ply is assumed to be 0.127mm. The total thickness of the patch is 2.032mm and 1.016mm, 
respectively. 

The move limit is set to be 5.0 degree in each new iteration step. The move limit will be 
reduced to half of former value if the new objective stresses after changing the fiber 
orientation are higher than that in the last step. The calculation will be terminated if the 
difference of objective stresses between the current step and last step is smaller than 10~ 4 or 
the move lim it is smaller than 10" 1 . 

Numerical results show that the optimum values of the adhesive stresses and optimal lay-ups 
are dependent of the initial values of fiber orientations. This clearly indicates the complex 
dependency of the objective adhesive stress on the ply lay-ups. To obtain a global optimal 
solution to the problem, evolutionary types of methods, such as Genetic Algorithm (GA), 
may be necessary. However, it is extremely computational expensive for the present problem 
because structural symmetry can no longer be used in each finite element analysis. 

To illustrate the algorithm, a very simple example is considered, in which all the fiber are 
aligned in one direction, i.e., uni-directional composite patch. In this case, there is only one 
design variable, i.e., the angle between the fiber direction and the longitudinal direction of the 
shell. 

Figure 5-1 to Figure 5-4 depicts the maximum adhesive stresses with different stacking angle 
and thickness of uni-direction patch of circular or square shape. From these figures, it is seen 
that the fiber direction affects the maximum adhesive stress in a very complicate manner, 
especially for the square patch. It is noted that the maximum peel stress has several stationary 
values ranging as the fiber angle from 0° to 90° for both shapes of patch. The maximum Von 
Mises stress has several stationary values too for the square patch. Table 5-1 gives the 
optimum maximum adhesive stress and optimum lay-ups of composite patch with the 
maximum peel stress and Von Mises stress being the objective function, respectively. 
Although 0° fiber can lead to a minimum peak value of Von Mises stress, the peel stress is 
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very high. Once again, it is necessary to determine the proper lay-ups according to the failure 
mode of adhesive layer. For example, if peel stress causes failure, the optimum fiber 
orientation is 50° for the circular patch, 60° for the square patch of 2.032mm thick and 0° for 
the square patch 1.016mm thick. 



0 to 20 30 40 50 60 70 80 90 

Ply-angle (degree) 

Figure 5-1 The maximum adhesive stresses with different stacking angle of uni-direction 

circle patch (thickness=2.032mm) 



0 to 20 30 40 50 60 70 80 90 

Ply-angle (degree) 

Figure 5-2 The maximum adhesive stresses with different stacking angle of uni-direction 

circle patch (thickness=1.016mm) 
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0 10 20 30 40 50 60 70 80 90 

Ply-angle (degree) 

Figure 5-3 The maximum adhesive stresses with different stacking angle of uni-direction 

square patch (thickness=2.032mm) 



0 10 20 30 40 50 60 70 80 90 

Ply-angle (degree) 

Figure 5-4 The maximum adhesive stresses with different stacking angle of uni-direction 

square patch (thickness=l.016mm) 
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Table 5-1 The optimum maximum adhesive stress and optimum lay-ups of uni-direction 
composite patch with different objective function 


Shape of 

patch 

Objective 

function 

Lay-ups 

Peel stress 

(MPa) 

Shear stress 

(MPa) 

Von Mises stress 

(MPa) 


Initial* 

[901 6 ] 

5.04 

18.53 

18.76 

Circle 

Peel stress 

[501 6 ] 

2.87 

12.13 

12.17 


Von Mises stress 

[01 6 ] 

10.83 

8.05 

10.20 


Initial 

[90s] 

10.73 

17.64 

18.80 

Circle 

Peel stress 

[50 8 ] 

4.83 

11.08 

11.43 


Von Mises stress 

[0 8 ] 

7.38 

5.96 

7.37 


Initial 

[90 16 ] 

9.49 

22.79 

23.44 

Square 

Peel stress 

[601 6 ] 

9.32 

21.53 

22.19 


Von Mises stress 

[0i 6 ] 

10.54 

8.31 

10.30 


Initial 

[90s] 

12.48 

21.00 

22.20 

Square 

Peel stress 

[0 8 ] 

6.57 

6.32 

7.20 


Von Mises stress 

[0 8 ] 

6.57 

6.32 

7.20 


^Initial lay-ups of optimization 


6. CONCLUSIONS 

In this report, a parametric optimization analysis has been conducted for optimal design of 
composite bonded patch repair to curved surface. Sequential Linear Programming (SLP) 
method and one-dimensional direct search method are employed as the optimization 
algorithm in conjunction with a fully implemented mesh generation algorithm into which new 
features were incorporated. Three types of design parameter are considered independently, 
which are: (a) geometrical shape of bonded patch, (b) ply drop-offs contour of bonded 
composite patch, (c) fiber orientation of each ply in bonded composite patch. The present 
study demonstrates the following salient points: 
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(1) The geometrical shape of bonded patch to a curved surface affects the maximum 
adhesive stresses significantly. It is useful and practical to select an optimal shape of 
bonded patch to minimize the maximum adhesive stresses, which makes the adhesive 
layer more reliable and sustain higher loading. 

(2) It is assumed that there are only three out-of-plane stress components in the adhesive 
layer, i.e., peel stress and two shear stresses. These two shear stresses can be 
transformed into a single one shear stress. According the optimization analysis, the 
peel stress can be reduced much more than shear stress does, which can be seen for all 
three kinds of optimization analysis. This is because the applied load is mainly 
transferred by shear stress in the adhesive layer from the parent structure to the 
bonded patch. 

(3) Three kinds of objective function are considered, which are the maximum peel stress, 
the maximum shear stress and the maximum Von Mises stress. It can be seen that the 
optimal patches may be similar or vary when different objective functions are used. 
But it should be pointed that it is necessary to choose a proper objective function 
according to its actual failure mode to determine the final optimal patch. 

(4) An automatic mesh generation algorithm is fully implemented with some new 
features being incorporated into. The present mesh scheme can fully fulfill the 
requirement of stress analysis and optimization evolution for bonded patch design. 


(5) Nodal stress and Gaussian Stress may lead to similar and slightly different optimum 
patch when used as objective function being used. Due to the nature of stress 
concentration along the boundary of the patch, it is believed that accurate nodal stress 
should be used as objective function. 

(6) The introduction of ply drop-off in a bonded composite patch can dramatically reduce 
the peel stress. However, the Von Mises stress may also increase because of ply drop- 
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offs. Therefore, one should be careful to apply this method to reduce the adhesive 
stress. 

(7) The effect of the residual thickness of patch with ply drop-offs may vary with the 
stacking sequence and the shape of composite patch. 

(8) The effect of fiber orientation on the maximum adhesive stresses is very complicate. 
According to the analysis, there exist several peak values in the whole design range of 
fiber direction, even for a uni-direction composite patch. The current method, SLP, 
can only present a local optimal solution efficiently. However, because of the required 
computational time, it is unacceptable to search for a global optimal solution by 
available global optimization method. Consequently, only a local optimal solution is 
given in this report. 
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